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Abstract 



Vortex dynamics in superfluids is investigated in the framework of the non- 
hnear Schrodinger equation. The natural motion of the vortex is of cyclotron 
type, whose frequency is found to be on the order of phonon velocity divided 
by the coherence length, and may be heavily damped due to phonon radia- 
tion. Trapping foreign particles into the vortex core can reduce the cyclotron 
frequency and make the cyclotron motion underdamped. The density fluc- 
tuations can follow the vortex motion adiabatically within the phonon wave 
length at the cyclotron frequency, which results in a further downward renor- 
malization of the cyclotron frequency. We have also discussed applications on 
the dynamics of vortices in superconducting films. 



PACS numbers: 64.40.Mj, 67.40.Db, 67.40.Rp, 67.40.Vs, 74.60. Ge 
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I. INTRODUCTION 



Besides collective excitations of density fluctuations, superfluid systems can also have 
topological excitations called vortices. Vortices are associated with many important phe- 
nomena, such as quantization of circulation Kosterlitz-Thouless phase transition 0, 
mutual friction [^,0, and flux creep Although the static properties of vortices are 

relatively well understood, the dynamical side is still wide open to investigations. 

It has been a common practice to treat a vortex in two dimensions as a particle, and to 
describe its motion by Newton's law in the following form 

Myito = -2TThpro X z - r^ro, (1) 

where tq denotes the vortex position, and p is the 2D superfluid number density. The 
terms on the right hand side represent Magnus and frictional forces respectively, and My 
stands for the vortex mass. The Magnus force term was motivated from the behavior of 
vortices in classical fluids, and has recently been related to the Berry phase of the many- 
body wave function Its existence has been experimentally established for superfluid 
^He, but this has been a controversial subject for magnetic fluxes in superconductors 0. 
The friction term is considered to be a result of the interaction of the vortex with collective 
excitations, or with the core excitations in the case of superconductors [0,0. However, very 
little experimental information is available to test the theories on dissipation at low enough 
temperatures. The vortex mass has also been a topic of debate. For incompressible classical 
fluids, it has been customary to regard the vortex as massless. This point of view has also 
been adopted in some calculations for nucleation and motion of quantized vortices [llO[ . 



Another point of view is that quantized vortices should have flnite masses, roughly equal 
to the mass of the vortex core In references |[T^-[T7|, the vortex mass was found to be 



renormalized by the condensate motion to a value logarithmically divergent with the system 
size in neutral superfluids. For charged superfluids, the vortex mass turns out to be flnite, 
because the screening currents effectively replace the system size by the London penetration 
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length |[T^-[T7|. Some recent discussion on the nature of the vortex mass can be found in 



The size of the vortex mass may affect tunnehng and specific heat of a vortex lattice ||T8 |. 

In this work, we study the vortex dynamics based on the nonlinear Schrodinger equation, 
which has been used to model superfiuids in a semi-microscopic manner. This equation has 



been derived by Gross and Pitaevskii for a weakly interacting superfluid []T^^. In the 
appendix, we derive a nonlocal version of the nonlinear Schrodinger theory from Feynman's 
many-body trial wave function, which can take into account strong correlations in a super- 
fluid such as ^He. The derivation of the nonlinear Schrodinger equation for superconductors 



in the clean limit is given in Refs. ||2l|j24 . This equation contains both solutions for collective 
excitations of density fluctuations and for topological excitations of vortices. It provides a 
useful starting point for the study of vortex motion and its coupling to the collective excita- 
tions. An effective Lagrangian will be derived for the vortex coordinate and for the density 
and phase fields of the superfluid condensate. As we will see, the equation of motion for a 
vortex naturally contains the Magnus force, and the effect of the condensate motion. 

We will concentrate our attention on the cyclotron motion of the vortex. Such a motion 
is a natural solution of the phenomenological equation (1), with the cyclotron frequency 
given by 

if damping is ignored. We found that the cyclotron motion is also a natural solution of 
the equations of motion of the vortex based on the nonlinear Schrodinger equation. We 
found that the cyclotron motion is damped by phonon radiation from the vortex. For a bare 
vortex, the damping rate is found to be about the same as the cyclotron frequency. For a 
vortex with trapped particles in it, the damping can be much smaller than the cyclotron 
frequency, although the cyclotron frequency itself is also reduced. With the cyclotron fre- 
quency determined from the equations of motion for the vortex coupled to the condensate, 
we can use equation (2) to define the vortex mass. The frictional coefficient rj in equation 
(1) can be calculated from the rate of damping due to phonon radiation. 



We have carefully examined the adiabaticity assumption used in references [p^|- [T^ , [T6| Jri 
that the field of condensate phase follows rigidly with the moving vortex. We found that 
this is possible within a length scale of A ~ .^(1 + Me/Mc), where ^ is the coherence length, 
and Mc and Me are the masses of the vortex core and of the trapped particles. Outside this 
length scale, the condensate cannot follow the vortex motion. Because of this, we found 
that the logarithmic divergence of the vortex mass with the system size, which was found 
in references [p!^ - [l^ , [T6| , [T7| , is cut off by the scale of phonon wavelength,. 

The organization of the paper is as follows. In Sec. II, we will give the basic ingredients 
of the nonlinear Schrodinger theory, and obtain the static vortex solution. In Sec. Ill, we 
will introduce vortex motion, derive the effective Lagrangian and the dynamical equations 
for the vortex and condensate. In Sec. IV, we will solve for the condensate response to 
the cyclotron motion, and find the conditions when the adiabaticity assumption is valid. In 
Sec. V, we will study the low frequency motion of the vortex, which is possible when the 
mass of the trapped particles is large. In Sec. VI, we present our numerical results for finite 
frequency of the cyclotron motion. In Sec. VII we will present our conclusions. 



II. NONLINEAR SCHRODINGER LAGRANGIAN 



Our starting point is the nonlinear Schrodinger equation which may be derived from the 
following Lagrangian: 



L 



p 



(3) 



where m stands for the mass of the superfluid atom, V represents the interaction potential 
between the atoms, and p is the background superfluid number density. The sign of V is 
positive to represent a repulsive interaction. 

The natural length and time scales of the nonlinear Schrodinger equation are ^ = 
h/ [mVpY^'^ and r = h/(Vp). As will be seen later, ^ gives the length scale of a vortex, and 
r gives the time scale of a bare vortex in cyclotron motion. If we scale the Lagrangian by 
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the energy fi^p/m^ then the Madelung transformation il) = y^e*'^ puts the Lagrangian in 
the dimensionless form 



pS+lp\VSf + j-^\Vpf + l[p-lf 



£r. (4) 



The dynamical equations of the condensate follow from the variation of the action with 
respect to the phase and the density 

P + V-[pV5]=0, (5) 



The first equation is nothing but the continuity equation, whereas the second resembles the 
Euler equation of hydrodynamics of classical fluids. The normal modes of the linearized 
equations of the fluctuations around the uniform condensate describe the collective excita- 
tions of the system. The low frequency part of the spectrum is phonon like with velocity 
Cg = C,/t = ^Vp/m PJ20|] . However, the roton minimum is absent in the spectrum, because 
the short range atomic repulsion is not appropriately taken into account. In the Appendix, 
we give the derivation of an effective Lagrangian that may overcome this problem. 

The equations also allow vortex solutions. For simplicity, we consider a single vortex with 
unit circulation around the center tq in the x — y plane, with the phase S of the condensate 
wave function given by 

S = So = Q{r-ro), (7) 
where 9 = arctan[(?/ — yo)/{x — xq)] is the polar angle of r — Tq. The density po satisfies 

The asymptotic forms are easy to find: 



Po(r)= ' ' (9) 




A simple analytical expression that interpolate between the above asymptotic forms is given 

by a 

= TT2^r (10) 



III. MOVING VORTEX 

In this section we will derive the equations of motion for a moving vortex. We assume 
that the motion of the vortex has a small amplitude (not necessarily slow), such that the 
fields of density and phase of the condensate may be expanded as 

S = So{r~ro{t)) + S,{r,t), 

p = Po(r-ro(t))+pi(r,t), (11) 

where at any instant of time, 5*0 and po satisfy the static vortex equations (|^) and (|). 5*1 
and pi represent small corrections caused by the motion of the vortex. By substituting 
Eq.(pTD into the Lagrangian (^), we arrive at an effective Lagrangian for Si{r,t) and pi(r, t) 
as well as the vortex coordinates ro(t). There is no problem of redundancy in the dynamic 
variables; unlike the original phase field S{r,t), Si{r,t) is required to be single valued p4[ . 
Keeping up to second order terms in fq, 5*1 and pi, we find the new Lagrangian as 



L = - d'r 



PoSo + ^polV^op + ^[po - 1]^ + Po'S'i + piSo + poVSo ■ VSi + 
+^|V5olVi + [Po - l]Pi + PiSi + ^Po|V5ip + piV^o ■ V5i + ^pI 

' 1 V7 V7 I^POP , IVpoP 2 Pi V7 V7 , l^^ll' 



+ 7^NPo\ +1 — Vpo-Vpi ^r^Pi + o ^ Pi -T^Vpo-Vpi 



ipo 4po 8pt 8pg 4p^ 8po 

We will add an extra kinetic energy term of the form 



(12) 



^M^il (13) 

to simulate the situation with particles trapped inside the vortex core. It is understood that 
this may be an over simplification of the real interactions between the trapped particles and 



the vortex [^], but the main purpose of introducing this term is to provide a mechanism 
for controlhng the time scale of the vortex motion. The external mass Mg is measured in 
units of (m^^p), and a bare vortex is described by Me = 0. 

The resulting Lagrangian can be simplified considerably by using the equations for the 
static vortex, integration by parts and dropping the constant terms: 

L = - j d^r PoSq + pqSi + PiSq + piSi + ^polVS'i]^ + piVS'o ■ V Si + 

, IVpoP 2 Pi V7 V7 , ^ IV7 |2 , 1 2 

-Pi - T^Vpo ■ Vpi + — |Vpi| +-pi 



8p3 4p2 ■ 8po' ■ 2 

-fSipoVSo-hdi+^M.tl (14) 

The boundary conditions will be taken as pi = and VSi = as well as 5*1 = at infinity. 
Then, the line integral, which is taken around the boundary, only adds a constant to the 
Lagrangian and can be dropped. The dynamical equation of the vortex is then obtained by 
variation of the action with respect to tq. The equation, hnearized in fq, pi and Si, is 

- 27rro X z + / [Si^Po - PiV^o] d^r - MJq = 0, (15) 

where we have used the fact that Vro can be replaced by —V when it acts on (r — fq). The 
first term represents the well-known Magnus force. The second term shows the coupling 
between the vortex and the condensate, and the last term is the usual inertial force on the 
particles trapped at the core. Note that by comparing Eq.(|T5|) to Eq.(|l]), one can extract the 
phenomenological parameters, such as the vortex mass and coefficient of viscosity. Also note 
that, it is not possible to identify the vortex mass immediately from this equation without 
knowing the perturbations pi and Si. 

The linearized equations for pi and Si are 

Si + WSo ■ VSi + -^Vpo ■ Vpi - T^VVi - ^-^rfpi + ^Pi + Pi = i-o ■ V5o, (16) 
4po 4po 4pg 4p^ 

Pi + Vpo • VSi + poV^Si + VSo ■ Vpi = i-o ■ Vpo- (17) 

From these two equations, we notice that the dynamics of the condensate is driven by 
the motion of the vortex. In other words, the vortex is accompanied with backflow-like 



corrections, pi and 5*1, which on the other hand, act back on the vortex in a way determined 
via Eq.(|T5|). 

In the rest of the paper, we choose the origin of r to be instantaneously at tq. The 
dynamical equations (p!5|), ( |T6| ) and ( p!7D remain unchanged up to second order terms in ro, 
Pi or 5*1. Also, the boundary conditions are such that pi and Si are regular at the origin, 
while pi 0, VSi — at infinity (boundary). 



IV. THE EQUATIONS FOR CYCLOTRON MOTION 



We would like to consider solutions of the following form: 



ro = Re be-"^\Si + iy) 



(18) 



Pi = Re bF{r)e 



-i{ujt~e) 



(19) 



5*1 = Re ihG{r)e 



-i(ujt-e) 



(20) 



The constant parameters u and b, and the functions F{r) and G{r), are yet to be determined. 
The first expression shows that the motion of the vortex is of cyclotron type with size b and 
frequency u. The frequency is allowed to have an imaginary part to describe a damped 
cyclotron motion. A complex phase factor in b is immaterial, because it only affects the 
initial angle of fq. A simple angular harmonic analysis shows that the angular dependence 
of the condensate response has to be in the given form. In a general sense, this motion is 
the analogue of the massive branch of the helical vortex waves in classical fluids p6| , p7| . For 
simplicity the expression "Re" will be dropped in the rest of the discussion. 

Using the following expressions for the right hand sides of equations (|TB|) and (^) 



ro 



ro ■ fp'o = -ibujpQe 



-i(u)t-e) 



(21) 



we find that Eqs.([l5|),(p!q) and ([T7|) imply the following equations for F, and G\ 



^ r°° M 

- [rp',G + F]dr + -^u = l, 

2 Jo 2,71 



(22) 



4pg 4po 



F' F' 
F" + 



P'o _ Po+p'o/r 
4p|] 4p2 









1 


F + 


'uj ^ ' 









G = -, 
r 



(23) 



p'qG' + po 



F = -ojp'q, 



(24) 



where a prime denotes differentiation with respect to r. Note that the parameter h has been 
factored out from the above equations, and it is only to be fixed by the initial velocity of 
the vortex. 



V. CONDENSATE RESPONSE TO THE CYCLOTRON MOTION 

In this section, we consider the response of the condensate density and phase to the vortex 
motion. Analytic and semiquantitative results may be obtained by studying the limits of 
large (r ^ 1) and small distances (r <^ 1). The large r limit carries the information on how 
the system size may affect the dynamics, and the small r limit gives the contribution of the 
core where most of the variation in po occurs. These two limits cover the essential features 
of the dynamics. 

When r <^ 1, we can neglect uo and replace po by 2r^ (cf.Eq.(^) on the left hand side 
of Eqs. (|2^) and (0). After a straightforward but tedious calculation one can show that the 
solutions are of the following form 

F(r) = 4Aa;r^ (25) 

G{r) = Bur. (26) 

It is not difficult to find A + B = —1, although the individual values of A and B are 
undetermined yet. Our knowledge of their sum will be enough for the purpose of evaluating 
the integrand in Eq. (p2D . 
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When r ^ 1, we may use the approximation of po for large r as given in Eq.(^). Also, 
the first three terms on the left hand side of Eq.(^) is dominated by F, yielding 

1 



F + 



00 — 



G 



00 

r 



(27) 



Substituting this relation into Eq.(|21), we arrive at the following equation for G 



G" + -G' 
r 



(28) 



It is easily verified that a special solution is G = ^. To obtain the general solution, we note 
that the corresponding homogeneous equation for G is that of the Bessel functions of order 



V = vT+~2u;. We only keep the outgoing wave, which represents the radiated phonons, then 
the general solution is given by 



G = - + GHl^\oor), 



(29) 



where C is a constant and if^^^-* is the first kind Hankel's function. This implies via Eq.(p7D 



(30) 



A similar radiation process is also found for an oscillating object in a classical fluid |2B . 

The constants B and C must be determined by matching the solutions in the intermediate 
region r ^ 1. For semiquantitative purposes, we may regard the above solutions for small 
and large r to be valid up to the coherent length r = 1 (where po = 2/3). Then the continuity 
of G and its derivative yield 

2 



C 



ouHi'^\ou) - H^'\oo) 



(31) 



and 



B 



looH('y{oo) + H('\ 



^ uoHi^'^' iuj) - hI^\uo) 



(32) 



where a prime denotes differentiation with respect to the argument. 
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VI. LOW FREQUENCY SOLUTIONS 



It will be instructive to consider the case of small tu, which, as we will see later, may 
be achieved in the limit of large mass of the trapped particles. Using the expansion of the 
Hankel function for the limit of small argument, we can find the constants as C = —i'^uj, 
and -B = |. 

Depending on the assumption that the asymptotic forms of the solutions are valid in the 
whole regions above and below r = 1 respectively, the cyclotron frequency can be determined 
from Eq.(^). Obviously, this constitutes a rough estimate of the frequency, but we expect 
the qualitative behavior to be well reflected. Then, to lowest order in uj we have 



^ln(l/^) + ^? 



In the large Mg limit, we can solve for to as 



Me , . 

uj + -^u; = 1. (33) 



1 



This gives a vortex mass via Eq.(2) 

M, = M, + ln(^) + ^^. (35) 

Apart from the mass of the trapped particles, there is a hydrodynamic correction (the second 
term) which diverges logarithmically with Me. The imaginary part represents a decay rate 
of the cyclotron motion due to the radiation of phonons. It is seen that this decay becomes 
unimportant when Mg >> 7r^/2. 

To understand the logarithmic correction, we divide the r > 1 region further into two 
regimes: 1 < r < and r > In the former regime, we can use the asymptotic form 
of the Hankel's function, with the result: 

G = (36) 

r 

and 
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F = -. (37) 

r 

This gives a density response pi exactly the same as found using the adiabaticity assumption 
||r^- P^P!B| . We will therefore call this regime as the adiabatic regime. This regime has a size 
of one phonon wavelength A = ^, and can be very large when Me is large. The adiabatic 
following of the condensate with the vortex in this region is responsible for the logarithmic 
correction to the vortex mass. 

For r >> l/cj we can write the solutions as 

G=i-^.&(-+i), (38) 
r V 2r 

and 

F=l + ^c.,/^e^(-+f). (39) 
r3 V 2r ^ ' 

It is not difficult to show that to first order in tq, the condensate phase and density becomes 

^(r, t) = S'o(r) + 6y|^e*(-+f )e-*-*, (40) 

and 

p(r, t) = po(r) + z6^y|^e*("'-+5)e-''^*. (41) 

The first terms on the right hand sides are the phase and density of a static vortex at the 
origin, and the other terms are oscillations in the condensate that represent the phonons 
radiated by the motion of the vortex. Therefore, at distances larger than A the motion of 
the condensate is not an adiabatic following of the vortex. 

We now calculate the coefficient of viscosity corresponding to the damping. It is given 

by 

(42, 

where Q is the quality factor, and we have restored the real units. Using the previous 
results, we find Q = M^/Mf,, where = ^^'^mp is roughly the mass of the superfluid 
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that is expelled from the vortex core. Using the value of the three dimensional density 
P3£i ~ lO^^m"^ for superfluid He^, the viscosity per unit length of a vortex line becomes 



^7 X lo- 
ci 



Me 



kg/ (m ■ s). 



(43) 



In comparison to the viscosity induced by the scattering of excitations in superfluid ^He at 
temperature 1 K 0], the numerical factor is on the same order, but the factor (Mc/Me) can 
reduce it further. 



VII. CYCLOTRON MOTION OF FINITE FREQUENCY 

In order to extend our conclusions to the small as well as to the bare vortex 
case, we adopt a semi-numerical method. We use the approximate solutions found in 
Eqs. ([25|) , (p^) , (p^) and (|30D to find a numerical solution of the cyclotron frequency from 
Eq. (p2D . Our results are contained in Fig.(|l]), where solid lines are the real and imaginary 
parts of the cyclotron frequency from the numerical calculation, and the dashed lines are 
from the approximate expression in Eq. (P3|) . 

The agreement between the approximate and numerical solution for large Me is quite 
satisfactory. The estimate. Me >> 7r^/2, for the radiative damping to be negligible is seen 
to be appropriate. In general we see that as Me decreases, the magnitudes of the imaginary 
and real parts increase, and for Me = they are both approximately given by the time scale 
of the nonlinear Schrodinger Lagrangian 

h 

uJr,r ~ = (44) 

Thus, the motion of a bare vortex is heavily damped. Also, we notice that the imaginary 
part is smooth over the whole range, but there is a peak in the real part near Me/27r ^ 0.81. 
This can either be a superficial defect caused by the inappropriate handling of the solutions 
at r ^ ^, or a real physical phenomenon, the nature of which can only be substantiated 
by a numerical solution of the complete dynamical equations. Such an effort is currently 
underway. 
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VIII. DISCUSSION AND CONCLUSIONS 



In summary, we have tried to give a better understanding of the dynamics of vortices in 
superfluid systems. We have derived an effective Lagrangian from a nonhnear Schrodinger 
Lagrangian, and obtained the dynamical equations for the vortex coupled with the conden- 
sate phase and density. 

We showed that the natural motion of the vortex turns out to be of cyclotron type just 
like the one predicted by the phenomenological vortex equation in (|T|). There are three quali- 
tatively different regimes for the condensate response: the core regime (r < the adiabatic 
regime < r < A), where A is the phonon wavelength at the cyclotron frequency, and the 
radiation regime (r > A). Combined with the numerical result, this wavelength is roughly 
given by A ~ ^(1 + Me/Mc). When the mass of the trapped particles is large compared to 
the core mass, the cyclotron frequency is low, and there is a logarithmic correction to the 
vortex mass due to the adiabatic following of the density and phase fluctuations in the large 



region of adiabatic regime [|T2|. Similar results have been obtained by Wexler and Thouless 
|29|, and Arovas and Freire using different methods. The phonon radiation damping is 
negligible here as was expected in an earlier work of Niu, Ao and Thouless . 

For a bare vortex, the cyclotron period is on the order of the time that a phonon travels 
a coherence length, and the vortex mass is on the order of the mass of the fluid that can 
occupy the core. The adiabatic regime is essentially empty, and the logarithmic correction is 
absent. However, due to the large density of states of the phonons at the enhanced frequency, 
radiation damping is heavy and may completely overshadow the cyclotron motion. 

In order to experimentally observe the vortex cyclotron motion, one has to create a 
situation of small damping. The present theory is too crude to exactly tell whether a bare 
vortex in superfluid ^He should be overdamped or underdamped, because the nonlinear 
Schrodinger equation does not treat the core structure accurately. One needs to carry out a 
microscopic calculation in order to pin down this issue. In the appendix, we briefly outline 
one such method. 
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Ion trapping in a vortex core can be a useful way of reducing the cyclotron frequency 
and the damping. We emulated this effect crudely by adding an external mass term to the 
Lagrangian. In reality, a trapped ion also expel the superfluid particles, making the core 
much bigger than the coherence length ^. For example, for negative ions (electrons) the 
expanded core size can be as large as 16A compared to a bare core size of lA in superfluid 
^He. In such cases the external mass in our theory should also include the expelled superfluid 
mass by the ions. The vortex can be driven either by an oscillating superflow, or by an 
ac electric fleld which acts on the ions. Resonances should be observed at the cyclotron 
frequency, which is on the order of where a is the radius of the hollow core. For the case 

^ ' ma'' 

of a negative ion, this gives u; ^ 6 GRz. 

When we consider charged superfluids, we must deal with the complication due to the 
existence of yet another length scale, the London penetration length. However, for thin 
superconductor fllms, the London penetration length is quite large, and vortices are very 
similar to their counterparts in ^He. For a bare vortex, the adiabatic length is about the 
same as the coherence length, so that there is no logarithmic correction due to the adiabatic 
following of the condensate. The cyclotron frequency should be given by the flux quantum 
times the background 2D charge density divided by the vortex mass: 



which only depends on the electron mass and the coherence length. Using the relation 
between the coherence length and the superconducting gap A, ^ = the cyclotron energy 
is found to be smaller than the gap by a small factor equal to A/e/, where e/ is the normal- 
electron Fermi energy. The cyclotron motion may be resonantly excited by an ac superflow 
as we discussed above for the case of superfluid ^He. For Al thin films, in the very clean 
limit, ^ = 1.6 /xm, we have u — MHz. The resonance should show itself as a peak in the 
resistance-frequency plot. 

The finite temperature effects can also be foreseen qualitatively along the lines of the 
present analysis. Primarily, since the core size increases with the temperature, the vortex 
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mass should become larger, and one expects to see the resonances at lower cyclotron fre- 
quencies. Another fact that contribute to this effect is the reduction of the Magnus force 
due to the reduction of the superfiuid density at finite temperature. 

Similar to the mechanism of phonon radiation damping in neutral superfluids, there is 
plasmon radiation damping in the superconductor case, but the effect is not strong because 
of low (zero for 3D) density of states of plasma at the cyclotron frequency. The cyclotron 
motion may also be damped by the excitations of normal electrons in the core. This effect 
might be strong due to the coincidence of the cyclotron energy with the excitation energies 
of normal electrons in the core, but the exact nature of this coupling has to be established 
in a microscopic calculation. 
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X. APPENDIX: A NONLOCAL NONLINEAR SCHRODINGER THEORY 



Assume that we are given the correct ground state |0) of a superfiuid in the absence 
of vortices and collective excitations. Then to a good approximation, the dynamics of the 



superfiuid may be described by the Feynman wave function ||3T| : 

exp (E,(^5(r,)+a(r,)))|0) 



1/2 



(46) 



(0|exp(2E,«(r,))|0) 

where 5* and a are time-dependent real functions, and the denominator is a normalization 
factor. The essential feature of the Feynman wave function lies in the factor of products of 
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single particle wavef unctions, whose position and time dependence may generate superflow, 
density fluctuations, and quantized vortices. 

The dynamical equations for 5* and a may be obtained variationally. The time-dependent 
Schrodinger equation can be regarded as the Euler-Lagrangian equations of the following 
Lagrangian 

L = ^(^|^|^)-(^|/f|^). (47) 

where H is the full microscopic Hamiltonian containing the kinetic as well as interaction 
energies of the superfiuid particles. If we substitute the Feynman wave function into the 
above expression, then, after some mathematical manipulations, the following Lagrangian 
may be obtained for the fields S and a: 



p(r) 

where the density p(r) is defined by 



^(r) + l|V5(r)|^ + l|Va(r)|^' 



dV, (48) 



^ (0|exp(2/a(rOp(rOrfV)p(r)|0) 
' (0|exp(2/a(r)p(r)d2r) |0) ^ ^ 

with p being the density operator. The last two equations provide a self-contained set. Incor- 
porated with the works on the ground state wavefunctions , this procedure constitutes a 



microscopic derivation of an effective Lagrangian that takes care of the short distance effects 
in a better way than the nonlinear Schrodinger Lagrangian. 

The usual nonlinear Schrodinger Lagrangian is recovered if we take |0) as the simple 
product state of zero momentum. In this case, the density p(r) is simply exp(2a(r)). In the 
general case, Eq.(^) still resembles the usual nonlinear Schrodinger Lagrangian by having 
identical terms involving the phase field S. The term involving a is generally a nonlocal and 
nonlinear functional of p, but is independent of 5*. A full scale study of the implementation 
and implications of the nonlocal nonlinear Schrodinger theory will be presented in separate 
publications. 
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FIGURES 

FIG. 1. Plot of the real (upper curves) and imaginary parts (lower curves) of the cyclotron 

frequency as a function of the external mass, calculated by using the analytical solutions (solid) 
and by using Eq.(33) (dashed). The frequency uj is in units of r~^, and the external mass Me.jl'K 
is in units of ^pm. 
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